Anti-phase locking in a two-dimensional Josephson junction array 
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Abstract 

We consider theoretically phase locking in a simple two-dimensional Joseph- 
son junction array consisting of two loops coupled via a joint line transverse 
to the bias current. Ring inductances are supposed to be small, and special 
emphasis is taken on the influence of external flux. Is is shown, that in the 
stable oscillation regime both cells oscillate with a phase shift equal to n (i.e. 
anti-phase). This result may explain the low radiation output obtained so far 
in two-dimensional Josephson junction arrays experimentally. 
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Josephson Junction arrays are currently under consideration as tunable microwave radi- 
ation sources.S'i. After clarifying the basic principles of (ID) arrayslHll, there is a growing 
interest in two-dimensional (2D) arrays since the pioneering work by Benz and BooiiH^. 
According to basic estimates^ radiation output generated by this type of arrays is expected 
to be much larger than that from ID arrays. However, experimental results point just in the 
opposite direction. While linear arrays delivered output powers up to 50/iW, the maximum 
output power reported for two-dimensional arrays is around lOOhwJiH^ and usually much 
smaller. While the general estimates referred to above are surely true, this should be caused 
by the fact, that actually very few junctions are locked in-phase. 

There can be several reasons responsible for poor radiation output. Besides technological 
problems this can as well result from the fact, that the basic mechanisms of phase locking 
in 2D arrays, despite some interesting results on several aspectsiSHS, have not yet been 
fully worked out theoretically. It is well known though, that there is no phase locking in 
unshunted 2D arrays in the absence of external flux. A theoretical study of the influence of 
flux with a "master-slave-mechanism" by Filatrella and WiesenfeldE! led to the conclusion 
that external flux can indeed lead to a certain phase locking; however the definite value of the 
phase difference could not be determined by their method, and stability was not considered 
at all. 

Here, we start with a very simple model 2D array, consisting of two loops coupled via 
a line transverse to the bias current (Fig. 1). Despite its simplicity this model is difficult 
enough to show the essential features of larger arrays: It is truly two-dimensional with a 
possible external flux entering the loops and a inductance in the transverse line, as it is 
typical in the nowadays favored hybrid arrays. Our propositions are as follows: (i) Both 
junctions are considered to be identical, (ii) Self-inductance is taken into account while 
mutual inductances are neglected, (iii) Junctions are overdamped with /3 ~ 0. (iv) There 
is no external load, (v) Instead of working within the framework of the widely used first 
harmonic approximation we exploit a phase slip technique which has proven successful in ID 
arrays before!!!. Its applicabiliy crucially depends on the proposition that the normalized 



ring inductance 

/ = 2nI c L/<t> (1) 

(Ic critical current of the junctions, L ring inductance, $o flux quantum) is sufnciantly small 
(/ « 1). 

Josephson junctions are described by the RSJ equations for the Josephson phases <pjk, 

4> jk + sin0 jfc = i jk ({j,k} = {1,2}), (2) 

where the dot denotes differentiation w.r.t. the normalized time variable 

s = ^R N I c t (3) 
n 

(Rn'- junction normal resistance; all currents are normalized to Ic)- Normalizing the exter- 
nal magnetic flux $ according to 

(f = 27T$/$ (4) 

we have to respect two flux quantization conditions, 

012 - 0ii - <f - U = 0, (5) 

022 - 021 - V 9 + H = 0- (6) 

In the following the transverse current playing a crucial role in the coupling will be denoted 
by i (cf. Fig. 1). In strong coupling problems of this type it has proven useful to introduce 
sum and difference variables according to0Ei 

S fc = ^(0fc2 + 0fcl), (7) 

A fc = -(0 fc2 -0fei)- (8) 
In addition, we introduce the circular currents 



i£ = (i*2 - i*i)/2. (9) 



With the help of these variables the problem can be reformulated as 

t k + sinSfcCos A fc = i Q , (10) 

Ajfc + sin AfcCosEfc = i° k , (11) 

Al + A 2 -^ = 0, (12) 

A x - A 2 - Z(z° - z?) = 0. (13) 

This indicates, that the voltage sums of both loops are driven by the bias current z'o > 1, while 
the circular currents drive voltage differences. Further, Eq. (|i~2l) is the flux quantization for 
the whole array, while Eq. shows that differences in the circular currents spread the flux 
differences of the loops. The transverse current z can be obtained from 

z = i°-z° = ~(A!-A 2 ). (14) 

According to Eq. (|i~3"D it is just the combination li which causes the coupling between the 
cells. 



The system (|I0|)-(|T3"|) is treated perturbatively assuming the ring inductance I to be 
sufficiently small. To lowest order, the flux quantization conditions gives (the second index 
indicating the order of evaluation) 

A M = <p/2, (15) 

i.e., junctions within both loops oscillate exactly in-phase. The Josephson oscillations itself 
can be evaluated from flTU|) as 

E *,o = 77 + 2 arctan — - — tan ( ° S l ), (16) 
2 iq + cos(y3/2) \ 2 J 

where we introduced the flux-dependent autonomous oscillation frequency 



Co = ^-cos^/2)). (17) 

This already completes the lowest order solution for our problem; Eqs. (|TI|) are not required 
for evaluating the Josephson phases within this order, but determine the circular currents 



i° k = sin(y9/2) cos S fcj0 (18) 

with 

v Co sin(C s - 5 k ) 

C0SS M = --— 7 J7y\ 77 "FT" 19 ) 

i + cos(y?/2) cos(Cos - o k ) 
To summarize, in lowest order the junctions within each cell oscillate in phase independently 
of the value of the external flux, while the relative oscillation phase between the cells remains 
undetermined. 

Changing to the next order I 1 we start again from the Josephson phase differences fll2"|) 
and (|T3]), inserting the lowest order result (|HJ) on the r.h.s. of Eq. (|H|). From the two alge- 
braic equations arising the correction terms A^i can be easily evaluated, and the Josephson 
phase differences of the two loops up to the first order in I are given by 

Ai = | + ^sin(^/2)(cosS 2i0 -cosS li0 ), (20) 
A 2 = ^ - ^sin(^/2)(cosS 2i o - cosE 1)0 ). (21) 

From this result, one can read off the transverse current 



sin(^/2)(cosS 2j0 -S li0 ) (22) 



with the basic harmonic 



, 4Co sin(y/2) ( S^ + SA . ( 6 2 - S 1 \ 

We point out, that although i is proportional to \jl this factor cancels out because of Ai — A 2 
being proportional to I itself. Accordingly, the amplitude of the transverse current is the 
same independently of the inductance I. 

The most remarkable property of this type of "internal shunt current" is its vanishing 
for ip = and growing with the external flux (p. One should notice, that this behavior is just 
opposite to that of an external shunt current, which usually turns out to be proportional to 
cos(</?/2). The absence of any transverse rf current for ip = is however obvious: In this 
case the array is completely symmetric. 



For evaluating the Josephson phase sums of the cells we exploit the method of "slowly 
varying phase" which has proven useful in the study of phase locking in one-dimensional 
arrays beforeffiS. According to this method corrections are put into the phases 5k, 

S k = 5 k {s), (24) 

which are supposed to change adiabatically only (in comparison to the rf Josephson os- 
cillations) in time. In addition, we will allow for the possibility that the joint oscillation 
frequency C be (slightly) different from the autonomous frequency Co- With these assump- 
tions the voltage sums can be written as 

* Co(C-fc) 



io + cos(<^/2) cos(£s — 5k) 
Inserting (^Sj) into flUf ) and neglecting higher orders in I after some algebra we arrive at 



(25) 



Co(C - Co - S\) = (Z/4)?(s) sin^ + (l/2)i i(s) sin^/2) cos(Cs - 5 1 ), (26) 
Co(C - Co - 4) = ~(//4)»(s) sintp - (l/2)i i(s) sin(^/2) cos(Cs - 5 2 ). (27) 

Here, all the interaction terms proportional to / arising on the l.h.s. of Eq. (pTDj ) were 
transferred to the r.h.s. In this way, the combination li plays a similar role as a synchronizing 
alternating external or shunt currentl'IHS 

To proceed, we average over one oscillation period, considering 8% as roughly constant 
over this time interval. It can be shown, that only the lowest harmonic (|23| ) of i contributes. 
Evaluation of the mean values results in the evolution equations 

Co(C - Co- <<*!>)= I , C °'° sin 2 (^/2) sin(< 6 2 >-<6 1 >), (28) 

A*o + so J 

Co(C - Co- < S 2 >) = -l ( C0t " : sin 2 (vV2) sin(< 5 2 > - < 5 1 >), (29) 

A 2 + SOJ 

where < 5k > denotes the one-period average over 5k- Subtraction gives the reduced equation 
for the phase difference 5 = 5i — 5 2 , 

<5>=l . %0 sin 2 (y/2)sin < 5 >, (30) 
*o + Co 
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having formally the same structure as the RSJ equation describing an unbiased autonomous 
junction. It admits two phase locking solutions, 

< S pl >= and < 5 pl >= tt, (31) 

describing in-phase resp. anti-phase oscillations of the cells. Investigation of the stability 
leads to the Liapunov coefficient 

A = l—^-r sin 2 (yV2) cos < 5 pl > . (32) 
l $ + Co 

As a result, only anti-phase oscillations are stable against small perturbations. By substi- 
tuting (|31 ) into fl28|) one easily recovers that the oscillation frequency remains equal to that 



of an autonomous junction, i.e. 

C = Co- (33) 



To summarize, the following picture arises: From earlier result st3£il we know, that the 
two junctions within each strongly coupled cell are generally (except for (p w tt) aligned 
in-phase. In addition, according to (|31]) both junctions from cell one oscillate anti-phase 
relative to those from cell 2. Synchronization of the cells in this state is provided by the 



alternating current (|22|) , flowing through the joint transverse connection. It is obvious, that 
such a state will be non-radiating. In addition, our findings justify earlier results on missing 
phase locking in the absence of external flux, which within our framework can be explained 
by the marginal stability observed in (|32|) for (p = 0. 

All results described in this paper are in complete agreement with corresponding numer- 
ical simulations performed in parallel. These simulations show, that the observed anti-phase 
locking is not bounded to the case of small inductances treated analytically here, but is a 
general feature of this type of array. If this remains true for larger arrays, which is under 
investigation now, this might well explain the low radiation output obtained with two- 
dimensional arrays up to now. In addition, investigations are on the way on the interplay 
with an external shunt current. 
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FIGURES 



FIG. 1. The two-dimensional Josephson junction circuit under investigation. 
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